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Uolng entirely nonproprietnry and nonsecret sources of information, a 
design program for an integrated Stirling cycle cryocooler has been written on 
on IBH-PC computer. The program is easy to use and shows the trends ond 
itemises the looses. The calculated results have been compared with some 
measured performance values. In its present form the program predicts somewhat 
optimistic performance. Tha program needs to be calibrated more with experi- 
mental measurements. As has been done before, adding a multiplier to vhe. 
friction factor can bring the calculated results In line with the limited test 
results so far available. The program is offered as a good framework on which 
to build a truly useful design program for all types of cryocoolers. 
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i. Introduction 

When I received a request from a client to 
produce a cryocooler design program for the 
IBM-PC, which both they and I had available, I 
assumed that such a program would be readily 
available from a government agency with an 
interest in this technology. It would simply 
be a matter of adopting it to the client's 
computer. We found, however, that access to 
all these programs is restricted, usually, to 
U. S. companies who already have a contract 
with a government agency to make cryocoolers. 
It is perfectly legitimate to have such re- 
strictions. However, there needs also to be 
available c program that can be used in schools 
and by companies who do not have an official 
need to know. 

Since 1 had written two design manuals for 
Stirling engines, and since many of the equa- 
tions that were U3ed in these manuals came from 
earlier publications showing how to design 
cryocoolers, I have undertaken the process of 
producing a cryocooler design program ana Logons 
to my numerous Stirling engine design programs. 

The previous literature will be cited. 
The method of analysis will be explained in 
general. The specific arrangement of a Stir- 
ling cycle cryocooler for which the program was 
written is given. A full list of input values 
is presented with typical input values. A 
sample of the full calculated output Is given 
and explained. Limited test results ore com- 


pared with calculated performance and dis- 
cussed. Application areas for thi3 . type or 
design program are suggested. Finally, conclu- 
sions are drawn about the utility of this type 
of design program. 

2. Literature Review 

A full literature review on cryocooler 
analysis is not attempted. Only those antece- 
dent publications that have a bearing on the 
computer program described in this paper will 
be given. In 1968 the author lead a group of 
engineers at the Donald W. Douglas Laborator- 
ies, in Richland HA In developing a Stirling 
engine for an artificial heart. He developed 
our own design mutaod. Gradually we became 
aware of other methods of analysis. We had 
simple methods and very complicated and time 
consuming methods. Later the author now, at the 
Joint Center for Graduate Study, University of 
Washington and sponsored by NASA-Lewls, wrote 
two design manuals for Stirling engines [1,2J. 
Also a long IECEC paper outlined in detail this 
design method [3). As a result of these publi- 
cations this method of design has been used 
widely among those who had no access to propri- 
etary information. 

There were a number of prior publications 
which were discovered in a literature search 
which had an important bearing in selecting the 
equations that were recommended in the design 



program which was presented. Crouthamel and 
Shelpuk [4] gave ail the design equations for a 
VM cycle cooler. Zimmerman and Longsworth { 5 1 
contributed to our understanding of shuttle 
heat loss as did Rios (6]. Leo [7] supplied 
the equation for appendix loss. Corrlng [8) 
supplied the equation for matrix conduction. 
Note that many of these publications are in the 
cryocooler literature. 

Since the design manuals have been written 
the author has continued, to evaluate and per- 
fect the isothermal analysis. A 1980 IECEC 
paper shoved that up to that time the isother- 
mal analysis was as good as any other available 
{9J. A further extension of this isothermal, 
program for large machines was done in which 
the effect of adiabatic hot and cold gaB spaces 
was included for both heat engines and heat 
pumps. This effect is particularly serious for 
engines or heat pumps that operate over a small 
temperature ratio. The 1983 IECEC paper showed 
that the improved isothermal analysis agreed 
with some engine test data over almost tb'. full 
operating range for both power output and effi- 
ciency to within plus or minus 10 percent with 
absolutely no adjustments { 10 ] . Finally the 
isothermal analysis has been extended to pre- 
dict the operation of free-piston machines and 
search for an optimum design 111]. 

3. Method of Analysis 

In the literature of Stirling engine de- 
sign, methods are classified as first, second, 
and third order. First order methods use var- 
iations of the Beale equation 112] to predict 
what the power can be, given the displacement, 
speed, pressure and temperatures. This method 
has been extended to large cryogenic refrigera- 
tors by Walker [13]. It does not apply to 
miniature cryocoolers. The method suggested in 
this paper is second order. It will be defined 
in more detail later. Second order methods are 
the simplest methods that take into account all 
the dimensions that influence performance. 
Third order methods use nodal analysis to simu- 
late the physical process very closely. They 
are expensive to run because they use a large 
amount of time on the very largest and fastest 
computers. They are useful for special stu- 
dies, but axe not practical to be used in 
search routines to find the best design. 

The basic assumption of the isothermal, 
rucond order analysis is that at each point in 
the cycle the pressure throughout the working 
gas volume is the same for each instant in 
time. In reality there are pressure differen- 
ces between the different parts of the working 
gas due to flow friction. In a well designed 
machine these differences are small in compari- 
sion to the pressure changes due to expansion 
and compression. Neglecting the pressure dif- 
ferences at this stage greatly simplifies and 
speeds the calculation with little loss in 
eccuracy. 


It is also assumed that an effective gas 
temperature can be identified for each part of 
the working gas volume which wolds for the full 
cycle. Thi6 assumption is easy to defend in 
the regenerator where gas is in close thermal 
contact with the solid. It is not true and 
requires corrections { 10] for the hot and cold 
spaces of large machines. Careful measurements 
of the pressure and total gac volume during one 
cycle in a Stirling engine for artificial heart 
power showed that on thlo scale the isothermal 
assumption was better then the adiabatic as- 
sumption. The artificial heart engine is on 
the same order of slse as the miniature cryo- 
coolers. • ' 

Using the two above assumptions, it is a 
simple matter to calculate the pressure during 
the cycle usually using the perfect gas law. 
If the pressure of the working gas is plotted 
against the volume of gas in the cold space, a 
closed curve is produced. The area inside this 
curve is equal to the heat absorbed . by the 
working gas in one cycle. When multiplied by 
the frequency we call this the thermodynamic 
heat input. If the pressure of the working gas 
is plotted against the total working gas vol- 
ume, a closed curve is also produced. The area- 
inside tills curve is equal to the work required 
to be applied in a perfect machine during one 
cycle. When multiplied by the frequency, ue 
call this the thermodynamic power input. 

These basic thermodynamic values oust be 
modified by the many losses that occur in a 
real cryocooler. To the basic thermodynamic 
power must be added the flow losses of the 
working gas, chiefly in the regenerator matrix. 
Also the 3eal friction, the mechanical fric- 
tion, and the electric motor losses must be 
added to obtain an estimate of the pover re- 
quired. 

From the basic thermodynamic cooling ef- 
fect must be deducted all the thermal conduc- 
tion and radiation heat losses. Also the flow- 
losses in the cold part are converted to beat 
which must also be deducted. The matrix heat 
transfer is not perfect. Therefore, the gas 
returning back to the cold space is at a higher 
temperature than when it left and must be re- 
cooled. Finally, there is a loss involving the 
gap between the displacer and the cold finger- 
cylinder. Keat i6 transferred by this shuttle 
loss only when the displacer oscillates. 

Each one of these losses mentioned above 
is calculated by an equation. In the second 
order method there is assumed to be no interac- 
tion between the various loss mechanisms. In a 
miniature cryocooler accurate estimation of 
these loss mechanisms is especially important. 
We find that ia cost cases the mechanical and 
electrical losses are more important than the 
gas flow losses and the thermodynamic power in 
determining the electric power requirement. On 
the cooling side, we find that the thermal 
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5. -Input Values 


losses consume most of the thermodynamic cool- 
ing effect. Therefore the net cooling effect 
is. subject to relatively large errors because 
it is the small difference between large num- 
bers. 

In a miniature cryocooler there is so much 
heat transfer area for so little volume that 
one would be tempted to assume that the heat 
transfer is nearly perfect. If tests show that 
this Is so, calculation of the performance 
would be speeded and simplified. We did not 
make thi6 assumption. We calculated a heat 
transfer coefficient for the hot space and the 
cold space and also decided on a heat transfer 
area. By trial and error we found an effective 
cold gas temperature lower than the specified 
cold metal temperature such that the heat that 
can be transferred based upon the calculated 
area and heat transfer coefficient is equal to 
the heat that needs to be transferred. At the 
same time by trial and error we found on effec- 
tive warm space temperature that is higher than 
the specified warm metal temperature such that 
the heat that can be removed from the warm 
space is equal to the heat that must be re- 
moved. When the fractional error in the heat 
balance at both ends is less than the conver- 
gence criteria for two trials in a row, the 
calculation ends and the results are printed 
out. 

U. Cooler Geometry 

Although this program can be adapted to 
any type of Stirling cryocooler, it is particu- 
larly set up for an i: fegral cryocooler shown 
schematically in figure 1. An hermetically 
sealed electric motor operates two cranks. One 
crank operates the displacer in the normal 
manner. The other crank operates the displacer 
piston in the inverted manner. The cold finger 
containing the displacer is in the same plane 
as the compressor piston. . The twp cranks are 
offset 90 degrees to provide the proper phase 
relationship. The program takes into account 
the slightly non-sinusiodal nature of the pis- 
ton and displacer motion due to the crank geo- 
metry. 


Figure 1. Schematic of Cooler Geometry 



The program stores all the Input data on a 
disk file. The user has the option of starting 
with the set of input data which was used last 
which la on this disk file, or I’stng base case 
input data which is written into the program. 
The operator makes this decision to start the 
program. Then the main menu is displayed on 
the screen (table 1). The operator can start 
calculating the next case by keying in a 1 or 
look at one of seven submenus by keying in 2 
thru 8. 

Table 1. Main menu for input. 

BUSIEST 

1 — Sian csiocistioi cast ossa. 

2 — Prcrrsj osatroi taramns. 

3 — 6«iET QSS'Ati C3 CESbtletS. 

t — fiestas era cosasttioa rc£j ass tress ess a. 

5 — thSD 41® CI8ES35KJ. 

. 6 — Eretssreaor oresmica. 

7 — Csltj eioa keeks*. 

8 — tfcst crai. 063!. cm Mice. 

Y<rcr cfcaics is — <fer fcsffisr sad aaml -5 • «.- 


Program control parameters. Table 2 shows 
the menu for all the program control parame- 
ters. On the left hand side is the menu num- 
ber. In the aiddle is the definition. On the 
right hand aide is the current value. Mote 
that at the bottom of the table a Ec*mu number 
is called for. Tie operator puts in an integer 
between 1 and 6 (Ho decimal point.). If a real 
number with a decimal point, or an integer 
outside of the specified range is entered the 
program displays the appropriate errer message 
and asks for the menu number again. If 1 is 
entered the program displays the master, menu 
agein. 

Table 2. Subnenu for Program Control Parasenters 


Provna control rerissters: 

: 1 — n«oro to esia tsaa. 

2 — Coo versKxs ctitetia. — 

3 — Cos* nsstor csf i«d b7 ojsrstor. 3S 

A — Graphic ostia — 8=ao. 1=?«. 1 

5 — PrintiM ostio? — £>=30. l^vcs. P 

6 — febsr of ties sices s?r circle. ; 24 

7 — Pressure droo cnion r— — 1 


0 o Catcvloie ft® pressure droo corretiti*. 

1 = Cdlcourte frto pressure drop test. 

8 — febcfft loss option — 8 

8 = Caicoiite froo best tresefer correlstiea. 

1 = Calculate fro Pressure dr® test 
Mips fe?>»td3 aasioer. 

9 — fbltiPlier for friction factor 1.® 

fsnu mater? 1 

If 2 is entered, the program then asks 
that a real number be entered with s. decimal 
point. This convergence criteria in compared 
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with the fractional change from one trial to 
the next of the thermodynamic poser input end 
the thermodynamic heat input. This criteria 
must be a small positive real number. The 
program checks to see if the number input is 
between 0.0 and 0.1. If it is, the convergence 
criteria is changed to the new value, and the 
menu is displayed again with the new value in 
place. If the number that is input is not 
within the specified range, an error message is 
displayed, and the menu is redisplayed without 
change. 

If 3 is entered, the program asks for an 
integer. If this case number is negative, the 
program rejects it and asks for another. 
Otherwise the case number can be set to any 
number desired. The case number is automatical- 
ly indexed without operator attention. 

If 4 is entered, the grephle option can be 
changed. Some IBM-PC*s have the ability to 
draw graphs on the screen. Others do not. The 
program can be used either way. If the compu- 
ter does not have the ability to do graphics 
and the graphics option is set at 1, the compu- 
ter will probably hang up. The program will 
have to be restarted. 

If 5 is entered, the printing option can 
be changed. Since some computers are equipped 
with printers and others are not, the program 
can be used either way. If a printer Is called 
for but is not connected, the program will 
stop. 

If 6 is entered, the number of time steps 
per cycle can be changed. More time cteps 
leads to higher accuracy at the expense of 
speed of computation. Twenty-four has been 
found to be a good compromise. A formula has 
been developed to compensate for a small number 
of tine steps (l). However, is has not been 
used in this particular program. 

If 7 is entered, the option on the way 
flow losses are calculated can be changed. For 
the 0 option, the flow loss is determined based 
upon a correlation by Kays and London cs sim- 
plified by Martini [2j. For the i option, the 
flow resistance test used to qualify the dis- 
placer for use is employed instead. Other 
options can be added as needed. 

If 8 is entered, the option on the way the 
heat transfer coefficient for the regenerator 
is calculated can be changed. For the 0 op- 
tion, a simplified Kays and London correlation 
is used [2j. For the 1 option, the heat trans- 
fer coefficient 1 g derived from the flow fric- 
tion test using Reynolds analogy. Other op- 
tions can be added as needed. 

If 9 is entered, a multiplier on the fric- 
tion factor can be changed. If- the Kays and 
London correlation is used, this factor is 
multiplied by the calculated friction factor to 


give the friction factor used to calculate flow 
loss. If the pressure drop test data is used, 
the computed pressure drop is multiplied by 
this factor. 

Engine operating conditions. Table 3 
shows the menu for the engine operating condi- 
tions. The menu is shown twice to illustrate 
how it is changed. The inputs are self expla- 
natory. The program adjusts the working gas 
inventory for the first 4 cycles to make the 
time averaged working gas pressure equal to the 
input value. After that leakage continues to 
operate between and crank case and the working 
gas. An effort has been made to ask for the 
input values In customary units. They are then 
converted to a consistent set of units in the 
program. 

Table 3. Submenu for Engine Operating Conditions 


Essie? oasratiw conditions: 

1 — Rstora to tfiio esui. 

2 — fiver is? worsitw as pressure, esia Sfifl.ES 

3 — fetor ssssdi rra 12C3.P3 

4 — fetal tea»f store in told ssace. K Tft.65 

5 — fiettl temperature in esro sasce. f SSI© 


Hrui nosfeor? 4 

Vitus? (with decimal oolml61 
Ecsiae ooeratira conditions: 

1 — Return to sain i. 

2 — fiverese mors ins pressure, csio 4ES.0? 

3 — fetor ssrod. rta 1752-?? 

4 — fetal ueseratore in cold swee. K 0.63 

5 — fetal tco»ratore in *ssra e sate, f 

Menu tvw&sr? 


Pistons -- connecting rodo — crankcase 
Table 4 shows the menu which allows values 
connected with the mechanical motion to be 
changed. The same organization is used a6 with 
previous tables. There is enough space to 
clearly define each input value and give tire 
units. At present the volume of the crank case 
space is not known very trell and is not used in 
the calculation. Strictly, b never, it should 
be used in leakage calculi ' ' ^ns and in kinetic 
simulations of the cooler. 

Table 4. Submenu to change inputs related to 
pistons, connecting rods and crankcase 


Pi ’.tens and ccwsectinu rods csd crati case. 

1 — Return to tain cssnu. 

2 — Creak ansle. dssreos S3.® 

3 — Dicseter of cceprcss or piston, inch ,5£S 

4 — Seal diaoner at disf®t?r. inch ':87a 

5 — Raisa of compressor Piston cr»K. inch ,^33 

o — fed: us of displacer drive crastt, inch ,gSj=t 

7 — Unsth of compressor piston ccon. rod. inch .dies 

8 — lerath of diSPlntc-r drive cons. tod. inch l.Tlfi? 

S — Vo l pe? of ctcns car? space, co. inches „ 


Mno lumber? 1 
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Usro flow passages Table 5 shows the 
submenu which allows values connected with die 
warm flow passages to be changed. Most of the 
entries ore self explanatory by reference to 
figure 1. The wans displacer duct Is the hole 
thru the seal and bearing part of the displacer 
to allow gas to pass from the connecting tube 
to the warm end of the regenerator matrix. The 
number of velocity heads are used to account 
for the flow resistances of entrances, exits 
and bends In the gas flow passages. Standard 
handbook, values were used (13j. The dead vol- 
ume Is an estimate of the dead volume In the 
warm space In addition to that In the piston 
end clearance, connecting tube, and the warm 
displacer duct. 

Table 3. Submenu for Marta Flow Passage Values 


Uiro fleet oasts* SJ 

1 — Gityra to turn ejtw. 

2 — DtcjsMr of CBKESctiss wfc** icO> .1623 

3 — Pisssier of ears diselecer duct. icch .6® 

a — Effective cessressor sistoa cloaraace, inch .G3S 

5 — Usstti of coorttiBS tube. lech 2.2y$ 

6 — Uasth of earn tucaiscet dttt, inch .463 

7 — tictSxr of wiccifr t tssia <a ceaaxtios! 2.25 

8 — (laasr of velocity bs&fs in diePlacsr etect 2.24 

9 — Bs&i voles? fit taro eas of dies later, co. inch _ .6233 
(teas ostrter? 1 


Regenerator properties Table 6 shows the 
submenu which allows values connected with the 
regenerator to be changed. The length of the 
regenerator matrix Includes space for the two 
coarse screens on each end. The porosity of 
the matrix Is computed from the total volume of 
the fine matrix, the weight of the fine 
screens, and the density of stainless steel 
(7.817 g/cc) 

Table 6. Submenu for Regenerator Properties 


ftaaserstor ercrartieuJ 

1 — {Saturn to esio e*su. 

2 — Insid* Cirartcr of disoiecet criicdsr. inch .SK3 

3 — DiesTtcr of r®KM»rs«s>r estria. inch . 1263 

4 — Hire diestar is rescseretor estrii, Inch .C£3 

3 — Sauare ssa h of screens in rso. ostris. uirwfia— S23.E2 

6 — tfcahn of fin* sere? as io res. estri*. CS9 

7 — Letwth ot ressaerstor esrtria. itth 2.23S) 

8 — Ifeitht of firs ressnercWT acmrs. sroa 1. 3123 

9 — taster ot cesres end Bcrs*s* in sotri# 4 

10 - Thictewa of each coare? end screes, inch .0143 

11 - Test floa rats of 12 thro disal., std cc/ein. 7SQ.E3 

12 - Pressore arcs at fc-st fira rata Mi 1S.S3 

few) nsntor? 1 


Cold flow passages Table 7 shows the 
submenu which allows values connected with the 
cold flow passages to be changed. In this case 
there io a single orifice which squirts gas 
Into the cold space. The correlation by Hauser 
(16 | was used to calculate the heat transfer 
coefficient in the cold space. 


Table 7. Submenu for Cold Plou Passages 


COtd flOd Pd 5 S* 3 *$! 

1 — Ertiira to nsia eras. 

2 — Disaster of cold orifice, lech - — — . .{5753 

3 — Effect!*? disoificar ciearoce fit old cc& luos_ .G32 

4 — Uwcth of cotd ofifico, inch — — _ .Ci'i'3 

3 — KrsSer of velocity beads is cold cfifice 1.9 

ttfS» Boaer? I 


Heat conduction, seal, and motor proper- 
ties Table S shows the submenu which allows 
values connected with heat conduction, seal, 
and motor properties to be changed. All the 
Input values In this list, except for the Inch 
dimensions are difficult to determine. The 
emmissivitleo can vary over a wide range depen- 
ding upon the exact properties of the 6urfocos, 
However the heat loss thru the vacuum insula- 
tion Is small In any case, eo these errora are 
unimportant. The efficiency of the drive motor 
and the seal and mechanical frlcltlos are very 
Important. They must be measured by separate 
tests. In thl6 case the efficiency of the 
motor was estimated from the motor specifica- 
tions. The seal and mechanical friction uas 
determined by adjusting it so that the electric 
power demand ot the very beginning of a cool 
down test was about right. The apeed of the 
motor during this test wao not recorded and had 
to be estimated. In reality the 3eol and mech- 
anical friction should depend upon engine apeed 
and working gas pressure. More realistic fric- 
tion con be added when the information Is 
available. 


Table 8. Submenu for Heat Conduction, Seal end 
Motor Values. 


Heat cssdt'ttiCB. ejai, sad tartor vaisss* 

1 — totsira to cam esss. 

2 — Biflsjtor of Vaasa issalCTiw dsfiaer. inch .6J 

3 — Efficieetv of driw estor. i E&S3 

4 — Eaiesivity of com fiaosr .C*S9 

3 — Eflissivity of tststa ieslstita chstfctr -G33 

6 — Cc2?icsd cast toasts, «.>ackfcsc/9Si - , .£2gj 

7 — Saal erf £*ch23icsi fr ictica. ia«s 7.G25 

8 — Thiesros of cold fiesr tan, inch .CS3 


9 — Uts> betesra cold f lecsr eai I asJ dissiecor, ieoi .£223 
{few basket? 1 


6. Calculated Outputs 

During solution two outputs are possible. 
If the graphic option is on, a display Is shown 
on the screen like that shown in figure 2. 
This display is very usefui to show at a glance 
whether the solution is converging and vhere 
some of the problens In design are. The ellip- 
ses at the right show the gas pressure plotted 
against the total volume. The first cycle is on 
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top. The other 6 cycle are below as the gas 
presoure is adjusted. The ellipses at. the left 
are the cold gas volume plotted against the 
working gas pressure. The sine curve at the 
top Is the displacer motion for one cycle. The 
sine curve at the bottom Is the compressor 
motion for one cycle. The nearly horizontal 
line at the top left records the trials of the 
effective warm space temperature. Full vertical 
scale Is 0-400 K. Hote that most of the 
adjustment is made after the first cycle. The 
nearly horizontal line at the bottom left re- 
cords the trials for the effective cold space 
temperature. Note that there was little 
change. 


Figure 2. Graphical Display during Solution 



If the graphic option is eff. Table 9 Is 
displayed on the screen. The numbers in the 
first two columns are the numbers that are 
compared with the convergence criteria which is 
shown in the top line of Table 9. These frac- 
tional changes for cycle 6 and 7 were all under 
0.001. Note that the work in and the heat out 
are now only changing In the fourth significant 
figure. The last column, the gas inventory, 
has unusual units. It Is the gas inventory In 
gram moles multiplied by the universal gas 
constant, 8.314 J/(gmol*K). 

Table 9. Progress to convergence table 

Ccs*x«aco criteria is* .531® {er6a8 SI 


Cms feasts 

Omss 


tones 

EM 

Ties to 

tteib Fosm la 

Cool la 

In 

la 

Pressars 8«* Iavsot. 

Fret. 

Fra. 

its I9S . 

JoylffS 

fi?a 

fee c. m 

1 1.63KS9 

• SSS575 

.1S3ES7 

.654T& 

2.ES3 

tSSlS .B132S3 

2 .672583 

.1SS73 

.172221 

.S2371 

2.SK3 

LS315 .612333 

3 .5552747 

.ScS 53 

•1E3234 

.CH2S8 

zem 

1.353 .6120738 

4 .SSM2S 

.g5379 

.1S23S0 

.6385 1 

2 . sas 

1.S313 .6128512 

3 .EE31 

.231517 

.IBS® 

.S5523 

ztm 

1.5315 .8123333 

6 .B25333 

.(ESI 37 

.1S2743 

.C5S331 

2.6273 

1.S31S .6126377 

7 .GS333 

.633443 

.lessi 

.(5S5S7 

2.533 

1.3313 .812373 


Speed of Solution On the IBH-PC equipped 
with the 8087 coprocessor, and running Micro- 
aoft FORT?JiN-77, the solution time for the 
above same case with 7 iterations is 6 seconds 
without graphics and 20 seconds with graphics. 
This does not Include the time for printing out 
or displaying the solution. The display of the 
output to the screen takes 2.4 seconds addi- 
tional. . The printout takes much longer and 
depends upon the type of printer that is being 
used. Ulth some type of spooler, printing can 
lag behind. Meanwhile the operator may inves- 
tigate new possibilities using the Information 
available at the monitor. 

Calculated Performance At the end of the 
soLutlon, the table of calculated performance 
is displayed on the monitor and Is printed out 
as well, if the print option is on. Table 10 
Is the calculated performance for the Input 
values given in Tables 2 to 8. The first thing 
to notice Is that some losses are very impor- 
tant and some sre negligible. However, It is 
dangerous to neglect the negligible losses . In 
the design calculation. If you do your search 
for a better design will lead you into on area 
where these negligible losses ere no longer 
negligible. 

Table 10. Calculated Performance 

Kama! Eej. fesisvic r» Dtiriiss Cvcit Cmucoicr 
lateral tesstrr. «€S.0iU5© FESti3S53So* 

K587 !SKflfS®T. tflTTa CO&ICl 3FSET. !STf3 
Thsjsswsmatesr &51ei Thsr.BmCtailES iff. S.6177 


Cid. Otic. F.l— .fjfl Cia. tV7t. Frit. feat .6S2S 

Ess. Flea 1(23 _ J.23S5 Gks 8. Fric, Ksss-— .5133 

Dies. fracJ.F.U— .@55 7 Shmio (test las: .1257 
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On the power requirement side, cote the 
mechanical friction and the electrical aotor 
losses are the two big ones accounting together 
for 74.4 X of the power requirement. The flow 
loss thru the regenerator accounts for 5,5 % of 
the power requirement. All the other flow 
losses account for only 0.2 7., leaving 19.9 X 
for the necessary thermodynamic power. Obvi- 
ously the piece to start ulth improvements ic 
to find some way to reduce the mechanical and 
electrical looses. 

On the cooling effect side, it appears 
that only 28.2 X of the thermodynamic cooling 
effect survives the gountlet of losses to be- 
come net cooling effect. In this analysis it 
is assumed that all the cold orifice flow loso 
and half the regenerator flow loss smst be 
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deducted as heat loss from the thermodynamic 
coo Hog effect. 

The shuttle heat loss end the appendix 
loss must be considered together, because as 
one Increases the other decreases. Shuttle 
heat loss occurs by redial conduction across 
the gap between the cold finger cylinder and 
the displacer. Because of the longitudinal 
temperature gradient along the cold finger 
cylinder and along the displacer, when the 
displacer la at the warn end of its stroke it 
is colder than the cold finger cylinder all the 
uey along, so it heats up. When the displacer 
is at the cold end of its stroke, it is warmer 
than the cold finger cylinder, so it cools 
down. Each cycle this process moves the heat 
that has been transferred one stroke length 
toward the cold end. This important heat loss 
con be minimized by using low conductivity 
cylinder walls such as plastic or glass or 
ceramic. The loss can be minimized by increas- 
ing the gap and increasing the length of the 
regenerator. 

Appendix loss occurs as the gas is packed 
into and unpocked from the appendix gap between 
the displacer and the cold finger wall. This 
gap is closed off at the warm end with a sli- 
ding seal and is open at the cold end to the 
cold space. As the pressure increases cold gas 
la packed into this space. As the pressure 
decreases not so cold gas flows back into the 
cold space. Increasing the thickness of the 
gap reduces shuttle loss but Increases the 
appendix loss. 

A number of equations have been published 
for shuttle heat loss and for appendix heat 
loss. As far as Is known there has been no 
experimental confirmation published of any of 
these equations. 

The reheat loss should, for a cryocooler, 
be called the recool loss. Both due to pres- 
sure changes and flow, a large amount of heat 
must be transferred each cycle In the regenera- 
tor. Ho matter how good the heat transfer, the 
gas always re-enters the cold space warmer than 
when it left. The reheat loss car be made 
8 mall at the expense of making the flow heating 
or the heat conduction loss or both too large. 
A careful balancing of dimensions is needed. 

The temperature swing loss is a correction 
to the reheat loss to compensate for the heat 
capacity of the regenerator being less than 
infinite. In this case it is negligible, but 
it can be important. 

The last 4 losses are heat conduction or 
radiation losses that go on whether the displa- 
cer moves or not. By far the largest is the 
cold finger wall since this must be thick and 
metal to hold back the pressure and contain the 
helium. The displacer' wail conduction is much 
less because this is made of plastic. The dis- 
placer matrix conduction is even smaller be- 


cause the metal screens are divided and make 
contact only in a few spots. Vacuum insulation 
loss la negligible. 

Record keeping The operator can continue 
to experiment with the computer program to see 
how different Inputs affect the output. He can 
do this without waiting for the printer on line 
at all. However with a spooler of some sort, 
the printer can keep track of all the inputs 
and all the outputs for further study. It taay 
be running many cases behind what the operator 
is considering. After the output of table 10, 
the full input values are printed out. These 
values are printed out in the same order as 
they are shown In tables 2 to 8. 

7. Results and Discussion 

Figure 3 shows how the computed perfor- 
mance uithout any adjuetmentof the calculation 
proceedure with test results. The net cooling 
effect is plotted against the cold finger tem- 
perature. Note thut if the flow resistance 
data is used for both the flow loss and the 
reheat that the net cooling 1 b too optimistic 
by 0.2 to 0.3 watts. If the Kays and London 
correlations are used, the net cooling is 
optimistic by 0.12 watts. From the dsta that 
oo far has been obtained. It appears that when 
the cryocooler motor speed was not measured 
that it was running at 1600 to 1700 rpm. Fig- 
ure 3 shows that the calculated and measured 
performance shows the some trend. 

For the flow rote that results in o 10 
psi pressure drop thru the displacer, the Keys 
and London correlation calculated 18.4 psi. 
Also a number of investigators have found that 
the flow resistance calculated for steady flow 
must be multiplied by a factor to make the 
measurements and the predictions for csiginea 
agree more closely. At one time the author used 
a factor of 2.9 to make the calculated output 
agree with measurements for both helium and 
hydrogen 117]. Tew [18] used a factor of 4 for 
hydrogen and 3 for helium to make his calcula- 
tions agree with test results. A recent test 
reported by Taylor and Aghili [14] shoved that 
oscillating flow in tubes increased the fric- 
tion factor by a factor of about 3. Possibly 
finding the right multiplier between reversing 
flow and steady flow resistance and the right 
way to apply this multiplier could bring the 
experimental and test results into agreement. 

Figure 4 shows the first attempt at adjus- 
ting the program to fit the current experimen- 
tal data. When the flou resistance test dsta 
are used, a mutipller on the computed pressure 
drop of 3-. 3 to 3.6 is needed to bring the 
computed results in line with the measurements. 
The measurements are consistant in that they 
show a trend with very little scatter. The 
computed results derived from flow- resistance 
test data do not show quite the same trend. 
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When the Kays and London flow correlation 
la used, the multiplier only has to be 1.25. 
Also the trends in the data and the adjusted 
calculation are the same. The correlation 
shows a transition from pressure drop propor- 
tional to flow rate to pressure drop propor- 
tional to the square the the flow rate in '-'a 
smooth manner over a factor of about 10 in flow 
rate. Possibly the pressure drop is not di- 
rectly proportional to the flow rate. 

These exercises in adjustment have been 
successful in showing that adjustments can be 
made. Houever, when we have more information, 
particularly the motor speed, for each point, 
our conclusions may be quite different than 
they are now. However, the isothermal method 
has an important advantage in that it itemizes 
the losses and therefore gives guidance about 
where to make improvements. 

8. Application Areas 

Although the specific cryocooler design 
program described in this paper is for an inte- 
gral Stirling cryocooler. The same programing 
concept can be applied to split Stirling and 
free-piston Stirling machines. A free-piston 
engine has been calculated using this type of 
analysis [11] . The speed for calculating each 
case is fast enough so that once the design 
method is calibrated, an optimization search 
can be programed -to search for the best design 
automatically. A thorough optimization search 
is already programed [11]. 

9. Conclusions 

Using only open sources of information, a 
cryocooler design program has been writ*--r. 
which gives reasonable results. It requires 
Inputs of. cooler dimensions, measured flow 
losses, mechanical friction, and electric motor 
losses. Since the current error between calcu- 
lated and measured performance is not large, 
and since the calculated performance shows the 
same trend as the measured performance, we 
expect that this computer program can be adjus- 
ted to model a real cryocooler quite exactly. 
The program needs to be adapted and customized 
for each user since there are so many ways to 
build a cryocooler. Also skill is needed in 
making the proper adjustments so that the pro- 
gram will be accurate over a wide range of 
design options. In itemizing the losses the 
program shows clearly where the greatest im- 
provements in design can be made. 
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Figure 3 Comparison of calculated cooling 
effect without adjustments with measurements. 



Figure 4. Comparison of calculated cooling 
effect with adjustments with measurements 
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